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The creation and destruction of entanglement between a pair of interacting two-level detectors 
accelerating about diametrically opposite points of a circular path is investigated. It is found that 
any non-zero acceleration has the effect of suppressing the vacuum entanglement and enhancing the 
acceleration radiation thereby reducing the entangling capacity of the detectors. Given that for large 
accelerations the acceleration radiation is the dominant effect, we investigate the evolution of a two 
detector system initially prepared in a Bell state using a perturbative mater equation and treating 
the vacuum fluctuations as an unobserved environment. A general function for the concurrence is 
obtained for stationary and symmetric worldlines in flatspace. The entanglement sudden death time 
is computed. 

PACS numbers: 03.65.Yz,03.67.-a,11.10.-z 



I. INTRODUCTION 



It has been shown fH [51 in algebraic quantum field theory that the vacuum state of a relativistic field is highly 
entangled. As by definition there are no real particles in a vacuum, entanglement defined here refers to the entangling 
effect two spacelike separated detectors would experience when placed into the vacuum. This inherently non-local effect 
is exhibited by the Feynman propagator which is non-vanishing between spacelike regions; a phenomena attributed 
to the exchange of virtual particles [3 - 5 . 

With an understanding of vacuum fields in inertial frames it is possible to progress to entanglement with the vacuum 
as observed by highly relativistic and non-inertial observers. Many such studies have now been performed [6Vll6j. 

In the case of linear acceleration one would expect the entanglement to be reduced by the interactions of the 
detectors with the Unruh bath [XT'. However for the accelerated two- level particle detectors shown in figure [l] it was 
found that 

\X\^e'^A (1) 

where a is the acceleration of the detectors and J7 the energy of the level separation. The requirement that an initially 
unentangled state become entangled can be expressed as \X\ > A. Thus, in the setup of figure [l] the detectors are 
always entangled. 

One can interpret A as the transition probability for the two-level detector to flip from ground to excited state, it 
does not contain any interaction terms between the detectors and therefore does not provide any entanglement. As 
the acceleration is increased, the transition probability per unit time approaches infinity. This is not surprising since 
the level flipping is enhanced by the Unruh radiation. What is surprising is that the amplitude \X\ (per unit time) 
also diverges. In other words, as the acceleration is increased, the cross detector vacuum interactions are enhanced. 
Since in this setup the detectors are separated by twice the inverse of the acceleration it is not clear whether it is the 
acceleration per se that increases |X| or if the increase is merely an artifact of the detectors being closer together. We 
are thus motivated to explore other systems of non-inertial motion to clarify this issue. 

The obvious place to start would be to enlarge the number of independent parameters in the linearly accelerating 
case, so that the distance between the detectors at the origin could be varied independently of the acceleration. 
However, this is not the simplest case; one flnds that the correlation function between detectors depends on both 
T — t' and r -|- r' in a non-factorizable way. Ideally we would like to find non-inertial motion such that the single 
detector correlations are stationary [19l |20] and the cross detector correlations can be factorized into functions of 
T — r' and r + t'. 
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When the non-inertial motion is not strictly Unear, the corresponding "Unruh" radiation seen by the detectors is 
known as acceleration radiation |21j . Typically on a stationary worldline the spectrum of vacuum fluctuations is not 
exactly Plankian but dependent on the motion of the detector [HI [23] . 



f 




FIG. 1: Interactions with a massless scalar field of particle detectors coUinearly accelerating in opposite directions. 
The amplitude X depends on field correlations that exist between the spactime events of the two detectors while the 
amplitude A depends on the field correlations that exist between events that occur at different times along the 
worldline of a single detector. At the origin the detectors are separated by a distance |. 



Regarding the physics, obvious experimental difficulties prevent the experimental interrogation of the vacuum in 
the linear case. Furthermore, the constant linear acceleration dilates the typical timescales involved to astronomically 
long magnitudes. However, for a particle following a circular path the acceleration 

is increased by a square of the Lorentz factor leading to a dramatic enhancement of the acceleration radiation effect, 
what is more since the speed is constant the time-dilation can be kept moderate; these points were first discussed 
in the context of measuring the Unruh temperature using single circulating electrons in [24] . a recent elaboration on 
the depolarisation of electrons in storage rings and it's connection to acceleration radiation can be found in [^ I26| . 
There is also the possibility of experimentation with electrons brought to high circular accelerations in a Penning trap 

For the reasons mentioned this paper will investigate the quantum information between two two-level particle 
detectors rotating diametrically opposite along a circular path with angular frequency uj, see figure[2| The symmetry of 
the detectors' motion is found to greatly simplify the problem. Importantly, both correlation functions are stationary. 
Furthermore, since the detectors have equal speeds the energy splittings and interaction timescales are not Lorentz 
shifted relative to each other, as would be the case if one of the detectors were kept at rest. The detectors at any 
given time are a fixed distance apart and have a uniform acceleration pointing toward each other. This provides an 
ideal setup in which to investigate the effect of acceleration on vacuum entanglement. 

Section [n] will investigate the significance of the vacuum entanglement by adiabatically switching on the interaction 
using a gaussian window function. Section [III| will focus on the density dynamics when the A interactions dominate 
starting with a maximally entangled initial state and investigating it's evolution and the time dependence of it's 
concurrence in the inertial frame defined by the center of motion of the particles. Throughout this paper the notation 
used is h = c = 1, Uj where j e x,y,z are the Pauli matrices and CTq is the identity matrix, and the metric is 
diag(-l-,-,-, -). 



II. ENTANGLEMENT WITH WINDOW FUNCTIONS 



The Hilbert space of the two circularly orbiting qubit particle detectors, A and B, and a massless scalar field, (j>, 
is written as H = C'j^ (E) Cg (E) H^. The Hamiltonian of the detectors, scalar field, and their interaction is given by 
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FIG. 2: Double helical structure of the particle detector worldlines, with radius R and angular frequency uj. 



H = Ho + + V, where, 



Ho = -^ri(cr^ (g) 1 + 1 ® 



V ^ ij{T){a^ ® 1 + 1 ® a^)<j){x{T)) (3) 

in which H^ is the free Hamiltonian for the field 0, and Q is the energy difference between the two states of a qubit 

particle detector [5H]. For convenience the interaction is turned on at time r = using the Gaussian, ?7(r) = 7706 ^ 
window function [101 El] , where ^ is a fixed parameter used to control the interaction time and 770 is the coupling. So 
far all of the above have been defined in the instantaneous rest frame (IRF) of the detectors, but it is more convenient 
to work in the inertial frame at the center of rotation and so we define^ , 

V'-"^, e=lC, a'=" (4) 
7 7 7 

and the worldlines for the two detectors are, 

XA — {t, R COS ivt, R sin uit,0) (5) 
xb = {t,—Rcosojt,—Rsmu}t,0) (6) 

where, 



R= and uj = (7) 



The positive frequency scalar field correlation functions, 



In fact a' defined here is not the coordinate acceleration in the inertial frame which would be given by ^ = Ruj^ = 7^a. Instead 
a' should be thought of as the change in coordinate velocity with respect to the change in proper time of a clock comoving with the 
detector, i.e., a' = ^ . 
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for the paths considered are, 



D+{xAit),XBit')) = - 
D+{xA{t),XAit')) = - 



1 



1 



(9) 



47r2 {t - t')^ - (2i?)2 cos2(a;(i - t')/2) ' 

1 1 
4^ (t - i')^ - (2i?)^ sin2(w(t - t')/2) ' 

and equivalent expressions for the correlation functions with A ^ B interchanged. 

Second order perturbation theory is used to calculate the final out density given it's initial state, p(— oo) — \ IDdi |. 
One begins with the Schrodinger equation for the density operator, p{t) — —i[V{t),p{t)]. Integrating over p{t), and 
substituting the result back into the right hand side, then integrating a second time gives, 

ft' 



p(oo) = p(-cx)) - i / [V{t'),p{-oo)]dt' 



dt' [V{t'),[V{t"),p{-oo)]]dt'\ 



(10) 

Assuming that the detectors and field are initially uncorrelated so that p{—(X)) = ps{—oo) (8" p^{—oo) and p^{—oo) = 
and solving ( 10 1 leads to the final out density. 



Ps(oo) = 



/ X 
OAF 
Q Y* A 

\X* 1 - 2A , 



where 



A 



X 



Y 



dt' / dt"f^'{t')ri'{t")e~'^''^''-'"^D+(xA{t'),XA{t")) 



dt' / dt"7^' {t')r^' {t")e'''''^''+'"^ D+ {xA{t'),XB(t")) 



2Re 



— oo 

t' 



dt' 



dt"rj'{t'W{t")e-^'''^''-'"'>D+{xA{t'),XB{t")). 



(11) 

(12) 
(13) 
(14) 



A necessary and sufficient condition for inseparability of (11), and also for entanglement, is that the negativity 
(defined as the absolute value of the sum of the negative eigenvalues of the partial transposed matrix) be greater than 
zero [5UH32| . As discussed in the introductory example, this amounts to |X| > A. 

Equations ( 12 1 and ( 13 ) can now be solved for the circular correlation functions ([9|. The amplitude A can be solved 
using contour integration, one finds poles distributed in the complex plane satisfying the equations 



zi — /3sinzi, 
Z2 = -/3sinz2. 

These equations are numerically solved in figure [s] An approximate solution for yo can be obtained^ 

yo 

Defining the dimensionless variables 

R 

and performing the contour integration we find 



y — fl^ a — a^. 




472 V r yoiyocothya - 1) 



y 



+ e-2yo!/VSerfc 



(15) 
(16) 

(17) 
(18) 

y 



E 



Zp{l — Zp cot Zp) 



-^""yy^erici-tzpj- +y) + e^''^''^^^ eTfc{~ 



'iz„\l y) 

a 



(19) 



^ Comparison with the exact numerical solution shows that this approximation is only valid in the rclativistic regime, we found a useful 
ansatz in the /3 — > limit was = ai {/3)arccosh(/3-l) where 1 < < \/3. 
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FIG. 3: Due to the symmetry of equations ( 15 ) the poles are symmetric about the real and imaginary axes. We 
define the independent poles to have positive imaginary part and negative real part, and label them in order of their 
distance away from the origin. All poles are first order except for zq which is a second order pole. 



where the gamma factor 

7^ = ra + 1 

has been kept explicit for presentational purposes. Also the poles are dependent on /3 which is 



13 = + 

ra 



-1/2 



(20) 



(21) 



in terms of the dimensionless parameters. 

The real part of X can also be solved using contour integration, one finds poles distributed in the complex plane 
satisfying the equations 



zi = /3 cos zi , 

Z2 — — 13 cos Z2- 

These equations are numerically solved in figure [4] An approximate solution for xq can be obtained 



(22) 
(23) 



(24) 



but the remaining poles must be determined numerically. The imaginary part of X can be computed directly (without 
the use of contour integration) adding both parts together we find that 



X 



1 7fo 



\ 



e ° (erfi(a;o^) - i) 
2xo{l + Xq tana;o) 



Im 



E 



zpeu{2<'=\4'='} 



e ^ {l + icr&{zpy^)) 
Zp{l + Zp tanzp) 



(25) 



By increasing the size of the window function, we can obtain the transition probability per unit time found in 



lim = — » . 



(26) 



^ The choice of hmit on the left hand side can be understood in the following way: the authors in 1241 define the transition rate per unit 
time by dividing the transition probability by d^(t + t') which corresponds in the case with a window function to dividing by a 



factor of limj_>oo ^ /+^ e ^ 
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FIG. 4: Due to the symmetry of equations ( 22 ) the poles can be divided into sets of four which are related to each 
other by complex conjugation and negation. We define the independent poles to have positive imaginary part and 
negative real part, and label them in order of their distance away from the origin. 



On the other hand we find that 

lim X = 0, (27) 

this can be understood by observing that the rate per unit time of X oscillates about zero and in the long time limit 
averages to zero. 

It is also satisfying that by replacing, a = and taking the the slow motion limit our equations reduce to 

the inertial case [TIJ)"^: 

2 

lim A = ^ (e-y^ - v^yerfc(y)) , (28) 
p^a 47r V / 

lim X _^le-'-^-^^ (crfi(r) - ^) . (29) 
We can now investigate the large acceleration limit while keeping the radius fixed, we find 

lim A = (30) 

a-i-oo Sa/StT 



hm \X\ = 0.42-^*^—— (31) 

(k) (k) 

where the decimal number in \X\ results from summing the contributions of the zl and poles up to k = 10 . 

Thus at fixed radius the cross detector correlations are reduced like the inverse square root of the acceleration. 
Furthermore, the A amplitude diverges proportionally with a. Both of these processes combined make it very difficult 
for any entanglement to be gained by detectors with large acceleration. The entanglement region is plotted in figure 

El 

Having shown that circular acceleration suppresses the production of vacuum entanglement, in the next section we 
will investigate the loss of entanglement due to the acceleration radiation. We calculate the entanglement lifetime of a 
pair of detectors that are initially entangled in a Bell state assuming for simplicity that the cross detector interactions 
responsible for the creation of vacuum entanglement can be ignored. 



* Upon making the correspondence in notation 2R — >■ L and ^ a. There is also the appearance of the imaginary term in X not present 
in | 10| . we have confirmed via private communication with the authors that the expression as we have it here is correct. 



7 




FIG. 5: Plot of the entanglement volume defined by the condition \X\ > A. This plot includes the pole contributions 
to equations (25 1 and (19) up to fc = 10. As the acceleration increases the entanglement volume is squeezed toward 

the a — y plane, thus at large acceleration the detectors only entangle for very small radii. In the opposite limit 
(a^ — >■ 0) the general features of the inertial case (TU] can be seen; a line in the r — y plane separates entangled from 

non-entangled regions. 



III. ENTANGLEMENT DYNAMICS 



In this section the dynamics of the two entangled detectors in the presence of a real scalar field are investigated. 
The method follows that given in [llj . Initially, the density operator of the two detectors are prepared in a maximally 
entangled state, called the Bell state pBeii = \{1 1 + ctx Ux — (Jy ay + az iz), which corresponds to the pure 
state (I ii) + I tt))/"\/2 written as a density matrix. As mentioned in the last section, since the vacuum entanglement 
is suppressed by acceleration it will be safe to ignore cross detector interactions. Each detector can then be considered 
in isolation before adding the two together to solve for the total density operator. For ease of notation we suppress the 
factors of i^ao so that for example cr+ should read cr+ ® (Tq- Utilising the Hamiltonian in equation ([s]), the equation 
of motion to second order for detector A in the Schrodinger picture is 

Ps = ^i[Ho,Ps]~Vo'{y] I dt'e'^'^''-'\a,G+ps-a+psO,)D+{xA{t),XA{t')) 

+ V / di'e-'f^'(*'-*)(a,a_p,-c7_p,a,)£'+(xA(t),a;A(t')) + ^-c.). 



(32) 



The correlation functions result from tracing over the unmeasured vacuum field, where Tr\p^(^(^'\ reduces to 
£'+(x^(t), a:A(i')) is the same as equation (|9). 

In order to simplify the second order equation (321) it is convenient to define /i = (/+ + /_)/2, /2 = — /_)/2, 

di = 1, and d2 = 0, in which 



/± = f dt'e^''''^''-*'>D+{xA{t),XA{t')). 
Jo 



(33) 



As this integral is sharply peaked and drops of quickly outside the domain of integration the Markovian approximation 
can be used to extend the integral to [0, oo]. By extending the domain of integration in this way the real part of the 



integral equation (33 1 is even, and can be solved by closing a loop integral in the upper (lower) half of the complex 
plane for (/_). The imaginary part is odd, diverging at the origin, and will need to be removed. 
After some simplification, the new notation gives the equation, 

Ps = -i[Ho,Ps] - Vo^ ^{(fjdt + fidj) {{a^aJ,ps} - 2ajPs<T,) + {fjd, - /*<ij)[cr,(Tj, p,.]|. (34) 

i,3 
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The last term, {fjdi — f*dj)[ai(7j, p^], when expanded is non-zero only for j = 2 and becomes — ilm(/+ — /_)[cr^, Ps]/2, 
which contains the divergent part of the integral (33). To remove the divergence the free Hamiltonian is redefined as 

/ 2 — 

H' — Hq + -^Im(/-|. — I-)az; the master equation now reading as, 

/ 2 

= -i[H', ps] - \ Y.^fjd, + f*d,)[{a,a,,p,) - 2G,p,ai) (35) 

where the energy separation between the states of the detector is renormalized as w'. 

Further simplification can be achieved by switching to the interaction picture in the basis a± to get. 

Pi =='7o^Re(/_)(2cr_p/o-+ - {(T+cr_,p7}) +77oRe(/+)(2(T+/?/(T_ - {cr_o-+,p/}) 

+ + r )fT+P/a+ + 77oe"''"'*(^- + lX)G^piG^. (36) 

The last two terms are rapidly oscillating due to the new renormalised energy, and can be dropped in the rotating 
wave approximation [33]. Following the same procedure for detector B and adding the two together results in the 
Lindblad form '3T of the total master equation, 

g=^(2i,pL]-{L]i„p}) (37) 
i 

where the operators are, 

Ly=^Q^Vic{IJ){G_®V) ^2 =?7o\/Re(7+)(CT+® 1) 



is = ??ov/Rc(/_)(l®(T_) L4 = r;ov/Rc(/+)(l®a+). (38) 

In order to solve the above system of differential equations the density operator is expanded via p = ^ TijUi ® Cj, 
and substituted back into the Lindblad equation. As the coefRcient, r^-, is symmetric when interchanging indices it 
is sufficient to look only at the top diagonal half of elements, which are, 

foo = f,, = 2?7^'(Re(/_ - /+)ro, - 2Re(/_ + /+)r,, + Re(/_ - /+)r,o) 

•tqx = -?7oVoa;Rc(/+ + /_) foz = 277o^(rooRe(J_ - /+) - rozRe(/+ -I- /_)) 

roy = -Vo^royRe{I+ + /_) r^^y = -2?7Q^Re(/_ + I+)r^j;y 

Txx = -2?yo^Re(/+ + I-)r^^ r^.^ ^ ryo^(2Re(/_ - /+)r^o - 3Re(/+ + /_)r^^) 



= -2ri'Q^ryyKe{I+ + I^) f^, = 77^,'(2Rc(/_ - /+)r2,o - 3rj,.Re(/+ + /_)). (39) 

Solving for r with the density operator initially being in the Bell state, the only non-zero elements are, 

''00 = ^ '''xx = ~Tyy ~ ^ ^ 

\5{l-e-"^-) r,, = i^2(l-2e-*/^^ + (l + l/<52)e-*/^i), (40) 



where 



= H'Re(/_+/+) (41) 
T^-i = Tf V2 (42) 



are the relaxation and dephasing times respectively. For convenience we have also defined 

^ Re(J_ - /+) 
Re(/_ +/+)■ 



(43) 



47r V 4^3 

Re(/+)=Re(/_)-^l)'. (44) 



The integrals (331 are solved by expanding the denominator in powers of {t' — t) to give, 

Re(^-) = 7Z + ' 
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FIG. 6: The diagonal elements of the density operator for a given acceleration, a, are plotted as a function of 



rescaled time t 



sit 



Substituting the results in (40) back into the expansion ps = ^ VijOi ® Oj gives, 







+ (l + ^(l_e-*/^^))2) 






4 






+ (l-<5(l-e-*/^^))2) 


Posit) 


= P3o(i) = 


2 



(45) 

The poo(i) component is the fraction of an ensemble of detector system pairs which would have both detectors in 
the ground state, while pssit) corresponds to both detectors being in the excited state. The remaining diagonal 
components, pii(i) and ^22(^)1 are the fraction of detector pairs with one detector in the ground state and one in the 
excited state. These matrix elements are plotted in figure [6] fo r various choices of acceleration. 
Taking the limit as time goes to infinity, from equation ( 45 1 the final out density reduces to, 



Psioo) = 2 n 1 



(1-6) 



1(1 + 5) 




(46) 



This can be expressed as the thermal state 



Ps[0O) 



Tre^'3'=ff^o ' 

implying the detectors have reached thermal equilibrium with the scalar field at an effective temperature given by 



(47) 



Toff — 

In ^ 

where the population of the ground state to excited state of an individual detector is, 



(48) 



P\ 



1-5^^ a 



(49) 
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the same as for a single detector in circular motion 
equations ( 41 ) and ( 43 ) as 



In terms of the effective temperature we can rewrite 



tanh 



(50) 
(51) 



The density ( 46 ) in the long time limit is clearly separable and thus all the entanglement initially shared by the 
detectors is eventually lost. To quantitatively understand how the entanglement evolves with time the concurrence, 
[35j . is calculated to give a value between zero and one, corresponding to no entanglement and complete entanglement 
respectively. Concurrence is given by. 



C{p) = max{Ai - A2 - A3 - A4, 0} 



(52) 



where {Ai, A2, A3.A4} are the ordered positive square roots of the eigenvalues of the matrix M — p {ay®ay)p* {ay®ay) 
From (45), the concurrence for the two detectors in circular motion is. 



C{p) = max{-(l - 5^){\ - e"*/^^) + ^(1 + 5^)t 



-t/Ti 



,0} 



(53) 



and is plotted as a function of both time and acceleration in figure [7j 



C[t',a/n] 




FIG. 7: The concurrence profile for the system of two detectors accelerating about helical worldlines, where again 
t' ^ rj^n/ 'J is the rescaled time and the acceleration is rescaled by the energy separation, 51. 

The curve along which the concurrence is zero is given by, 

= T2 In I I . (54) 

This behaviour is another example of what has been termed entanglement sudden death (ESD) [151 136]. It can be 
seen that for zero acceleration (6 = 1) the concurrence decays exponentially with time, but for finite accelerations 
there is a finite time before which the concurrence completely disappears. 

As the acceleration increases and becomes much larger than the energy difference between states, the time taken 
for concurrence to go to zero is 

iesd - 47ry31n ( — ^ ) (55) 
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The 7 factor is due to the time dilation in transforming from the instantaneous rest frame to the inertial frame. 
As mentioned in the introduction, this dilation is less sever than that which occurs in the linear case, where the 
transformation is essentially exponential [TT]. 

Careful consideration of the steps involved in arriving at equations ( 53 ) and ( 54 1 indicate that they will hold for any 
pair of two-level detector systems (assuming the cross interactions are negligible) moving on symmetric and stationary 
fT9l worldlines. 



IV. CONCLUSION 



This paper has considered two qubit detectors orbiting a common point in the presence of an unmeasured massless 
scalar field. This setup provides an ideal situation in which to study the vacuum entanglement for detectors undergoing 
uniform acceleration separated by some constant distance. 

We have found that acceleration suppresses entanglement between detectors by reducing the amplitude of the cross 



term like the inverse square root of the acceleration (31 1 and increasing the level-flip transition amplitude linearly 
with the acceleration (30). Although our results pertain specifically to circular motion we conjecture that the way 
in which these amplitudes depend on the acceleration is universal to all non-inertial motions of detectors interacting 
with a free massless scalar field in flatspace. This conjecture can be justified in the following way: at any given time 
the detectors are a fixed distance apart and have fixed acceleration vectors pointing inward, we expect the direction of 
these acceleration vectors to be irrelevant to the local distortions of the field and therefore to the entangling behaviour 
of the field on the detectors. 

Having established that the vacuum interactions could be ignored for large accelerations we then studied the 
dynamics of the two detectors initialized in a Bell state. For large acceleration compared to the energy gap of the 



detector levels the ESD time ( 54 ) was related to the relaxation time T2 ( 42 1 by 



tesd - 0.88ir2 (56) 

where T2 was found to be inversely proportional to the acceleration. This indicates that it might be easier to observe 
the radiation effect by measuring the entanglement directly rather than the components of density matrix. 

We expect our expressions for the concurrence to be useful in other studies of non-inertial entanglement when the 
worldlines are symmetric and stationary. In particular, these will hold for the case of linearly accelerating worldlines 
with an independent distance parameter ^. 
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